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Abstract. In this paper the response of the heliospheric termination shock to an 
incident interplanetary shock is examined. This paper is an extension of a recent study 
by Barnes (1993), which treated the analogous problem for an incident contact 
discontinuity. The termination shock is treated as a strong gasdynamic shock. The 
postinteraction configuration consists of a moving termination shock, a postshock 
contact discontinuity, and either a shock or rarefaction wave propagating the 
disturbance signal into the downstream medium. For a decrease in dynamic pressure a 
rarefaction wave propagates downstream, and the new termination shock propagates 
inward, while for an enhancement of dynamic pressure the termination shock moves 
outwards and a weak outer shock propagates into the downstream medium; speeds of 
motion of the termination shock are typically of the order of —100 km/s. The results 
are similar to those presented by Barnes (1993) indicating that the results of that paper 
are robust within the gasdynamic model, in the sense of being independent of the 
details of the initial disturbance. 


1. Introduction 

The equilibrium location of the heliospheric termination 
shock is variously estimated by various investigators [e.g., 
Suess , 1990; Baranov , 1990; Holier , 1989; Lee , 1988] but 
may confidently be placed in the range 50-200 AU, probably 
with significant variation over the solar cycle [e.g., Lazarus 
and Belcher , 1988; Lazarus and McNutt , 1990]. The pros- 
pect of a near-term encounter of various spacecraft presently 
in the distant heliosphere with the shock has led to increased 
interest in the properties and motion of the termination 
shock [ Barnes , 1991, 1993; Smith , 1991; Suess , 1993], We 
expect the shock to be in constant motion [ Barnes , 1993 
(hereinafter referred to as Paper 1); Belcher et al ., 1993; 
Suess, 1993]. 

In this paper we present a simple gasdynamic model of the 
motion of the termination shock in response to a change in 
the upstream solar wind conditions. This is an extension of 
Paper 1 , which analyzed the effect on the termination shock 
of an upstream jump in dynamic pressure due to a contact 
discontinuity (i.e., an increase (or decrease) in density with 
no change in the speed at the discontinuity). In that paper it 
was conjectured, but not demonstrated, that similar results 
would be obtained for initial disturbances more general than 
a tangential discontinuity. In this paper we explore this issue 
further in a model that assumes that the jump in dynamic 
pressure is manifested as either a forward or reverse inter- 
planetary shock. 

We construct a simple quantitative ideal gasdynamic 
model of the motion of the termination shock as in Paper 1 . 
The termination shock is assumed to be a strong shock; for 
simplicity we take it to be initially at rest with respect to the 
Sun, although this assumption is not necessary. Upstream of 
the shock is a discontinuous increase (decrease) in dynami- 
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cal pressure in the form of a forward (reverse) interplanetary 
shock, which eventually encounters the termination shock. 
The incident interplanetary shock and termination shock are 
both assumed to be planar and the motion is assumed to be 
one dimensional (in a Cartesian coordinate system) in order 
to make the calculation analytically tractable. While this 
approximation is inadequate as a global description, it is 
valid as a local and initial description of the interaction. In 
particular, it should give the correct near-term postinterac- 
tion velocity of the termination shock. Generalization to 
spherical geometry is desirable but will probably require 
numerical simulation beyond the scope of the present inves- 
tigation. The analysis shows that the postinteraction config- 
uration depends on whether the interplanetary shock is a 
forward or reverse shock. An encounter of a reverse inter- 
planetary shock with the termination shock results in an 
inwardly moving strong shock, identified as the new termi- 
nation shock, and a simple-wave rarefaction propagating 
downstream. In the case of an interaction between a forward 
shock and the termination shock the resulting configuration 
consists of two shocks with a contact discontinuity between 
them. The inner of these two shocks is an outwardly moving 
strong shock and is identified as the new termination shock, 
while the outer shock is a weak outwardly moving shock that 
carries the signal of the disturbance into the downstream 
medium. The inward and outward speeds of the new termi- 
nation shock depend on the magnitude of the change in the 
upstream dynamic pressure but are typically of the order of 
100 km/s. 

The results obtained in this paper are similar to those in 
Paper 1. This is to be expected, as the resulting configuration 
after the interaction of the termination shock with incident 
interplanetary shocks is the same as that after an interaction 
with an incident contact discontinuity. Apparently, the in- 
teraction is due mainly to the arrival of a jump in dynamic 
pressure and depends only weakly on the details of the 
variation. 
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idate this condition; in any case, in the formal calculations 
we admit all M w0 > 1. The upstream and downstream 
dynamical variables are related via the Rankine-Hugoniot 
jump conditions. 

We also suppose that initially there is an interplanetary 
shock Si somewhere upstream of S 0 , which eventually 
encounters S 0 . The interplanetary shock may be either a 
forward or reverse shock, that is, propagating either antisun- 
ward or sunward in the local plasma reference frame. Let 
M ui be the Mach number of S { as defined from the down- 
stream (antisunward) side of S f ; note that M ui is greater than 
(less than) unity if S ( is a forward (reverse) shock. Let the 
dynamical variables upstream of S; be indicated by the 
subscript 1 . Note that the density ratio pjpi is limited to the 
range (y - l)/(y + 1) < pjp x < (y + l)/(y - 1). The jump 
conditions allow us to express M ui in terms of this density 
ratio by 


*4 = • (i) 

(y+ I) —~(y - i) 

Pi 

However, since M ui and M u0 refer to the same plasma, 

( 2 ) 

where V ( is the velocity of 5/ as measured in A 0 . Combining 
these two expressions permits us to express V j!v u in terms 
of M u0 and the density ratio, that is, 


Figure 1 . Schematic representation in the x-t (space-time) 
plane of the interaction of an upstream interplanetary shock 
S i with the termination shock 5 0 . (a) If p x !p u > 1 after the 
interaction, there are two outwardly propagating shock 
waves S\ and S 2 with a contact discontinuity C x between 
them, (b) If p\lp u < 1 after thelnteraction, there is a single 
inwardly propagating shock wave, a contact discontinuity 
C { and a simple-wave rarefaction 2ft propagating into the 
downstream. 


The plan of the paper is as follows: in section 2 the 
assumptions of the model and details of the calculation are 
presented, numerical results are presented in section 3, and 
the conclusions are given in section 4. 


2. Formal Calculation 

The aim of the paper is to calculate the velocity ( V , ) of the 
new termination shock and the velocity (V 2 ) of the rarefac- 
tion wave or weak second shock that propagates down- 
stream after the interaction of an interplanetary shock with 
the termination shock. Consider a frame of reference A 0 in 
which the initial termination shock S 0 is static. We treat the 
plasma as an ideal gas whose ratio of specific heats is y. Let 
p w , v u » P w » and c u ~ {ypJPu) Xil i respectively, represent 
the density, velocity, pressure, and sound speed upstream of 
S 0 , and let p 5 , v s , p s , and c s represent the corresponding 
quantities downstream (see Figures la and lb). Let M u0 = 
vJ c u be the Mach number upstream of S 0 . We anticipate 
the M u0 » 1 , although a sufficiently dense population of 
interstellar pickup ions in the outer heliosphere might inval- 




M u0 


f (y + 1) — - (y - 1) 
Pi 


( 3 ) 


where the sign of the function sgn (x) = x/|x| distinguishes 
between the forward and reverse cases for shocks 5,. Using 
this expression in the velocity jump condition for permits 
us to give an explicit expression for the upstream velocity 

V t 9 


V \ Pu 


Pu\ V,- 


= Ml I = l 

Pi \ Pi/.*'* 

P_u_ 

Pi 


M 


u 0 


f (y + 1) — - (y - l) 

Pi 


( 4 ) 


Similarly, the square of the upstream sound speed C\ is 


(p u \ 

x Pi 


(y + l) (y - l) 

\p 1/ 

Pu 

Mlo 

Pi 


(y + 1) (y - 1) 


Pu 


( 5 ) 


Thus, for given v u and M uQ the character of the incident 
interplanetary shock 5, is completely determined by the 
density ratio pjp x . When 5,* encounters the initial termina- 
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Figure 2. Schematic representation of the postinteractior 
geometry after an interplanetary shock has encountered the 
termination shock. is the postinteraction terminatior 
shock and C x is a contact discontinuity, (a) \i p x !p u > 1 , th( 
transition between regions 3 and 4 is through a weak shod 
wave S 2 . (b) If Pi/p H < 1, the transition is through s 
rarefraction simple wave 01. 


tion shock 5 0 , the following configuration emerges as shown 
in Figure 2. There will be a new termination shock S x that 
will be in motion in the reference frame A 0 . Downstream of 
S 1 (region 2 in Figure 2) there will be a region of shocked 
plasma (characterized by dynamical variables p 2 , etc.) 
bounded by a contact discontinuity C x on the downstream 
side. Downstream of this discontinuity (region 3 in Figure 2) 
will be a region of material, originally shocked by S 0 , which 
has responded to the disturbance created by the collision of 
S 0 and S ( . The signal of this disturbance will be carried to 
the distant downstream region either by a second shock S 2 
(if S ; was a forward shock) or through a rarefaction simple 
wave 01 (if 5/ was a reverse shock). The far downstream 
region (region 4), which the disturbance signal has not yet 
reached, is characterized by p 4 = p s , v 4 = v s , p 4 = p s , and 
c 4 — c s . 

The detailed calculation of the response closely parallels 
the analogous discussion of Paper 1. It is convenient to 
perform the analysis in the reference frame A! in which the 
velocity in the far downstream plasma vanishes. The veloc- 
ity in this frame is denoted by u, where 


u = v — v 4 


( 6 ) 


and v 4 = v s is the velocity of the far downstream plasma as 
measured in the frame A 0 in which the initial termination 
shock Sq is static. 

The far downstream gas is characterized by the following 
equations, derived from the jump conditions: 


u 4 = 0 

(y + QA^o 

P4 (7 - 1 )M „ 0 + 2 P " 

_ Puvl [2yM* 0 - (y - 1)] 

P4 ~ y (y + l)A/„ 2 0 

The far downstream sound speed is given by 

= Av u 

where 

A = ^[ 2 yM„ 2 o-(r- D 3 [(r - 1)^0 + 2]. 

(y + \)M u q 

(9) 

The density and pressure jump conditions across the new 
termination shock 5 x are given by the standard jump condi- 
tions and are analogous to (17) and (18) in Paper I but with all 
the terms retained since we consider all M n > 0, where 
M\\ is the Mach number of region 1 with respect to 5, , and 
is given by 



(7) 


( 8 ) 


Af fi = 


(Ut-U,) 1 


( 10 ) 


while c i is the sound speed in region 1 and U ) the velocity 
of S ] measured in the frame Aj. 

Define the dimensionless parameter rj by 


Pi_ Pl_ _1 m Iq 

Pa Pu v\ 2yMl 0 - (y - 1 ) 


• [2y(«i - U { ) 2 - (y - l)c 2 ] (11) 

where we have substituted for M 2 , using (10). This expres- 
sion can be rearranged to yield 

U\-V x fp~u 2yA/ 2 o - (y - 1) (y - 1) c? 

v u VP i V 2 7 M lo + 2y v 2 u ^ 12 

where the sign is determined by the fact that region 1 is 
unshocked, so that u x > U j. 

Combining the expression for the velocity jump condition 
across 5 0 


y - 1 2 

yTT + ( 7 + i)m 2 0 


with (4), we get 


Vu V u 


V4 

Vu 


y + 1 
Pu 


1 - 


M 2 u0 


Pi 


- 1 


M 


«o 


(13) 


f(y+ 1) — -(y - 1) 

p 1 


( 14 ) 
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The velocity jump condition across S ] , which is analogous 
to (20) in Paper 1, provides the following equation after 
substituting ( 10 ) for M \ x : 

2 2 

7 ("1 - t/i ) 2 + («2 - «|)(«1 - l/i) - — 7 C? = o. 
y + 1 y + 1 


This yields the following relation: 
u i — U x y + 1 


Ml- "2 (mi-m 2 ) 2 16 /C|\ I 

■(—■ + V-^— ) 1,61 

where the positive sign must be chosen because u ] > U \ . 

Equating the two expressions (12) and (16) for U x lv u 
results in an expression for 17 in terms of M u0 , pj/p^ , y, and 
u 2 fv u , which may be expressed as 

\p~u y + 1 

J—vfi + c- — — 

V p 1 4 

" 7> 

where 


consider the transition between regions 3 and 4, which 
involves the second shock S 2 . The pressure jump condition 
together with definition of 77 (equation ( 11 )) yields an equa- 
tion analogous to (23) in Paper 1 , which can be rearranged to 
produce 

(y + 1)77 + y - 1 


where M 42 is the Mach number of region 4 with respect to 
the shock S 2 . For the configuration described above to 
persist w 3 < U 2i where U 2 is the velocity of shock S 2 , 
which means that region 3 will contain shocked gas with 
M 42 > 1 and 77 > 1 . Since u 4 = 0, M 42 - (U 2 lc 4 ) 2 , which 
yields the following expression for U 2 : 


V 2 y - 1 I y + 1 

— = A V~j — v + — 1 17 • 
V 2y v r ~ 1 


The continuity equation across regions (3) and (4) (with u 4 = 
0) yields u 2 = U 2 ( 1 - p 4 /pi)- Hence (21) and the jump 
condition for the density (see (26) in Paper 1) can be 
combined to give the following expression for u 2 : 


u 2 w 3 2 

v u v u Vy(y - 1 ) 


1 y + 1 

1 + T 71 

I y - 1 


2 yA^ 0 - (y - 1 ) 
2 ^ 


y - 1 /c 


2 y 

and Cj/l^ and are given by (5) and (14), respectively. 

The next task is to find an expression for u 2 as a function 
of 77 ; such an expression used in (17) will permit us to solve 
(numerically) for 77 , and hence complete the solution of the 
problem. The expression for u 2 will depend on whether 5/ is 
a forward or reverse shock. If S t is a reverse shock, a 
rarefaction simple wave $1 propagates into the downstream 
gas. has the properties outlined in Paper 1, resulting in the 
relation 


(y- 1) 


(c 3 - c 4 ) < 0. 


The expansion in 91 is isentropic (equation (42) in Paper 1), 
and substituting ( 8 ) for c 4 results in the following expression 
for u 2 : 


u 2 w 3 2 
V u V u (y - 1) 


(y- 1 )/ 2 -y 


where the dimensionless parameter A is given by (9). We 
now have an equation for u 2 (v u in terms of 77 and the known 
parameters, M uQi p x lp u and y. 

When Si is a forward shock (pi/p„ > 1), the configuration 
after the interaction of S { and S 0 will consist of a new 
termination shock S x , a contact discontinuity, and a second 
shock S 2 that propagates the changed upstream conditions 
to the downstream gas. Our analysis follows Paper 1 , and we 


Equation (22) is a function of 77 and the known parameters 
M u q , p 1 /p w , and y and can be used in (17) to solve for 77 . 

Note that allowable 77 must have an upper bound. Because 
region ( 1 ) is unshocked 


2 [ P u 

«2 < = — t y> — 

y + 1 \ pi 


1- 


— - 1 

y + 1 pi 


( (y + 1 ) (y - 1 ) 

Pi 


The term containing 77 in (22) increases monotonically for 
all 77 > 1. It then follows from (22) and (23) that 77 < 77 ^ , 
where 


V* ~ ^ 2 + 


1 + 1 + 


4(L 2 - 1) 
y + 1 

2 + V 

y - 1 


where 


\ 2 y(y - 1 ) 

L = ib. 

(y + UA 

In the limit M u0 — * oo we recover (29) of Paper 1 . 
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Table 1. Numerical Results for the Encounter of an Interplanetary Shock With the Termination Shock for M u0 = 100 


P\Pu 


CO 

II 





V 

Vi lv u 

V 2 lv u 

1? 

Vilv u 

V 2 lv u 

0.25 

0.452680 

-0.332201 

0.809170 





(0.445058) 

(-0.334253) 

(0.809017) 

(0.458497) 

(-0.354249) 

(1) 

0.333,... 

0.540191 

-0.262563 

0.809170 

0.554146 

-0.279340 

1.000125 


(0.532948) 

(-0.264454) 

(0.809017) 

(0.546483) 

(-0.280413) 

(1) 

0.50 

0.685817 

-0.164614 

0.809170 

0.697899 

-0.175103 

1.000125 


(0.679863) 

(-0.166073) 

(0.809017) 

(0.691624) 

(-0.176115) 

(1) 

0.666, . . . 

0.806006 

-0.095670 

0.809170 

0.814752 

-0.101713 

1.000125 


(0.801768) 

(-0.096655) 

(0.809017) 

(0.810287) 

(-0.102465) 

(1) 

0.80 

0.889668 

-0.052386 

0.809170 

0.895155 

-0.055666 

1.000125 


(0.887014) 

(-0.052981) 

(0.809017) 

(0.892352) 

(-0.056144) 

(1) 

1.00 

1 

0 

0.809170 

1 

0 

1.000125 


0) 

(0) 

(0.809017) 

0) 

(0) 

(1) 

1.25 

1.126380 

0.053068 

0.836753 

1.119508 

0.056021 

1.029364 


(1.122608) 

(0.052326) 

(0.835792) 

(1.115421) 

(0.055364) 

(1.028257) 

1.50 

1.237942 

0.096008 

0.860066 

1.224203 

0.101300 

1.054004 


(1.229910) 

(0.094495) 

(0.858258) 

(1.215303) 

(0.099888) 

(1.051812) 

2.00 

1.430263 

0.162994 

0.898291 

1.404497 

0.171971 

1.094569 


(1.411955) 

(0.159775) 

(0.894584) 

(1.382493) 

(0.168588) 

(1.089602) 

3.00 

1.746130 

0.257042 

0.956596 

00 

CO 

00 


(1.693384) 

(0.248693) 

(0.947027) 

(1.63578) 

(0.261583) 

(1.1435) 

4.00 

00 

00 

00 





(1.90866) 

(0.309229) 

(0.984619) 

(1.82564) 

(0.324419) 

(1.18166) 


V\iv u is the velocity of the postinteraction termination shock S x and V 2 /v u is the velocity of the second shock S 2 when p\/p u > 1, or 
the velocity of the head of the rarefaction wave when p\/p u < 1. The results from the calculations in Paper 1 are given in parentheses for 
corresponding density ratios. 


3. Numerical Results 

Equation (17) specifies 17 as a function of p\lp u where 
p\fp u is restricted to the range (y - 1 )/(y + 1) < p\lp u < 
(y + l)/(y - 1) and u 2 is given by (19) or (22) depending on 
whether the initial shock is a reverse or forward shock. The 
termination shock velocity U\ , follows from (12) or (16), 
with U 2 given by (21) for a forward incident shock and the 
speed of the rarefaction simple wave equal to c 4 for a 
reverse incident shock. A Galilean transformation to the 
frame A 0 , fixed with respect to the Sun, gives the velocities 
V x and V 2 of the shocks S x and S 2 respectively. For a 
reverse incident shock, V 2 is the velocity of the head of the 
rarefaction wave $1 and is given by c 4 + v 4 . 

In the limiting case 77 — » 0, which corresponds to the 
expansion of the gas of region 4 into a vacuum, p x , p 2 , p 2 , 
p 3 , and c 3 ail vanish, and the velocity u 2 reduces to 


u 2 w 3 2 

v u y - 1 


(25) 


In the strong shock limit (M m0 — > 00) with c x = 0 the limiting 
velocity of the new termination shock S x in the frame A 0 is 
given by 


v i = £ 1 + ^ = _ 

V U v u 



(26) 


The velocity of the head of the rarefaction wave is given by 
V 2 = c 4 + v 4 . Both V x and V 2 are functions of M u0 , y and 
p x lp u . In the limit 17 — > 1 , V x = 0 in the strong shock limit, 
and V 2 = c 4 + u 4 . In the limit 17 — > 17*, V x tv u = 1 and 
V 2 lv u = (3y - 1 )/(y + 1) for the strong shock case. 

First, we compare the results of the present calculation 
with those of Paper 1 (for which the incident disturbance is 
a contact discontinuity). In the latter case the density ratio 


p x ip u is permitted to range from 0 to «>, whereas in the 
present paper the permitted range of density ratios is 
( y - l)/(y + 1) < p\/p u < (y + l)/(y - 1). For this range 
of density ratios the results of the present calculation should 
be identical to those of Paper 1 in the limit M u0 — > 0°. This 
is confirmed by in Table 1, which gives results of the present 
calculation for M u0 = 100, and (in parentheses) the corre- 
sponding results from the calculation of Paper 1 . 

The results of the present calculation should differ from 
those of Paper 1 by a factor of the order of HM u0 , due to the 
appearance of a term of this order in (14). Table 1 shows that 
the calculated propagation velocities agree with those of 
Paper 1 to within a few percent, as expected. Also, the two 
values y (2 and 5/3) give qualitatively similar results (Figure 
3). Once V x is found other velocities of interest may be 
calculated: 


v 2 v 3 1 

V„ V u (y + 1 )M]\ 

• |[(y - l)A/u + 2] — + 2(M^ — 1) — i 

l Vu V U J 

c 2 ^[2yM?, ~(y~ l)][(y- 1)M?, + 2] /„, yA 

v„ (y + l)A/n \v u v u J 

(27) 

C3_ lyJXy + l) + (y - i)r?] *>2 

v u A y (y - 1) + (y + l)v Pu 

c -l = A-n (y ~ Wy — <1 

v u Pu 


where 
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Table 2. Numerical Results for Various M u0 With y = 5/3 


Pl/pu 


M u0 = 30 



o 

II 

o 

5! 



M u o = 1.5 


V 

Vi/w. 

V 2 lv u 

V 

Vllv u 

V 2 lv u 

V 

Vllvu 

V 2 /v u 

0.25 

0.470626 

-0.327169 

0.810719 

0.522813 

-0.310506 

0.824273 

0.902936 

0.080399 

1.398411 

0.33 

0.557185 

-0.257951 

0.810719 

0.606113 

-0.242864 

0.824273 

0.932947 

0.091849 

1.398411 

0.50 

0.699706 

-0.161068 

0.810719 

0.739057 

-0.149566 

0.824273 

0.968019 

0.089313 

1.398411 

0.67 

0.815841 

-0.093280 

0.810719 

0.843337 

-0.085537 

0.824273 

0.986073 

0.068956 

1.398411 

0.80 

0.895807 

-0.050942 

0.810719 

0.912797 

-0.046260 

0.824273 

0.994060 

0.045014 

1.398411 

1.00 

1 

0 

0.810719 

1 

0 

0.824273 

1 

0 

1.3984 11 

1.25 

1.135316 

0.054732 

0.840245 

1.161627 

0.059109 

0.859801 

1.364223 

0.107939 

1.509578 

1.50 

1.257081 

0.099393 

0.865603 

1.314206 

0.108212 

0.891517 

1.779623 

0.202937 

1.621942 

2.00 

1.474394 

0.170135 

0.908434 

1.609903 

0.188359 

0.948852 

2.847618 

0.372027 

1.866426 

2.50 

1.673846 

0.226233 

0.945425 

1.922765 

0.255320 

1.004740 

4.498870 

0.533206 

2.172024 

3.00 

1.877770 

0.275086 

0.981310 

2.315656 

0.318141 

1.069494 

7.590468 

0.711242 

2.624672 

3.50 

2.146327 

0.325049 

1.026040 

3.029991 

0.389863 

1.175600 

16.29955 

0.967821 

3.549101 

3.90 

2.903426 

0.400567 

1.140108 

6.308316 

0.504641 

1.555725 

81.53364 

1.631841 

7.176246 

3.99 

6.942233 

0.514081 

1.593834 

31.88829 

0.638652 

3.131357 

786.9706 

3.688561 

21.037996 

3.999 

35.30583 

0.629486 

3.236902 

252.1864 

0.820020 

8.311837 

7742.927 

9.953374 

64.734375 

4.00 

00 

00 

00 

00 

00 

00 

00 

00 

00 


V\!v u and V 2 fv u are the velocities of the postinteraction termination shock (Sj) and the second shock (S 2 ) if P\!p u > 1 or speed of the 
head of the rarefraction wave if p x < p u , normalized to the upstream velocity v u . 



and vj, Cj and r\ are given by (4), (5), and (11). 

Let us now consider lower values of M u0 . At distances of 
—50 AU the measured proton temperature is typically of the 
order of a few times 10 4 K, at least at low heiiographic 
latitude [Gazis et al ., 1994], and varies only slowly with 
heliocentric distance. No measurements of the electron 
temperature are available beyond —5 AU, but if we assume 
that the electron temperature is comparable to the proton 
temperature, then M u0 will be in the range —8-30 for a solar 
wind speed of 400 km/s. A population of interstellar pickup 
ions sufficiently hot and dense to amount to an appreciable 
fraction of the solar wind pressure would imply a lowering of 
this estimate. Table 2 shows results for M u0 = 30, 10, and 
1.5, with y = 5/3 (results for y = 2 are similar and not given 
here). The results for M u0 = 30 and 10 are similar to each 
other and to the results for M u0 = 100 (Table 1). The 
behavior of the new termination shock and the second shock 
or rarefaction wave does not vary qualitatively for values of 
M u0 down to —2.4. However, at values of M u0 lower than 
2.4 the new termination shock may move inward or outward 
depending on the value of the density ratio. When M u0 
reaches 1.6, the motion of the termination shock is outward 
for all permitted values of the density ratio, as shown in 
Table 2 for M u0 = 1.5. 

For larger M u0 values the inward shock speeds for density 
ratio decreases are typically larger than the outward shock 
speeds for comparable density ratio increases. However, as 
M u0 decreases, these speeds approach similar values until 
M u0 — 24. Then, for density ratios near unity the outward 
shock speed for a density increase is larger than the inward 
shock speed for a corresponding density decrease. The range 
of ratios for this occurrence becomes larger with further 
decreasing M u0 . For M u0 <2.4 the new termination shock 
S j can assume either inward or outward excursions depend- 
ing on the density ratio decrease as shown in Figure 4. In 


these cases the outward shock speed for p\fp u > 1 is always 
greater than the outward shock speeds for a comparable 
decrease in density ratio. The speed of the weak outward 
moving shock wave is significantly faster for lower values of 
M u0i as illustrated in Figure 4. 

A singularity arises in the solutions as p\/p u — > (y + 1)/ 
(y - 1). In this limit, (3H5) give infinite v\ f V h and C\ , 
corresponding to an infinitely strong incident interplanetary 
shock 5/. However, the speed V x of the new termination 
shock increases only very slowly in this limit; for example, 
even for the extreme case p\!p u = 3.99, for which the far 
upstream solar wind velocity v\ = 9.3 v u ( — 3700 km/s) the 
speed of the new termination shock is only V x = 0.82 v u 
(-330 km/s). 

Altogether the results of the present calculation are qual- 
itatively similar to the results of Paper 1, except for the 
extreme cases of small M u0 and effectively incident inter- 
planetary shocks of near-infinite strength. 

4. Conclusion 

The motion of the termination shock resulting from inter- 
action with an interplanetary shock has been studied. Our 
analysis shows that the postinteraction configuration de- 
pends upon whether the jump in density associated with the 
interplanetary shock, which is restricted to the range (y — 
l)/(y + 1) < p\lp u < (y + l)/(y - f l), is greater or less than 
1. In the case of a density increase (p, > p u ) the termination 
shock moves outwards, while for a decrement (pj < p u ) the 
resulting motion of the termination shock is inward for 
reasonable values of M w0 . For larger M u0 values (>24) the 
speed of the inward propagating shocks are slightly higher 
than those of the outward propagating shocks at comparable 
density ratios while for smaller M uQ (<24) the speed for the 
outward motion tends to be larger, particularly at density 
ratios close to 1. The speeds of both inward and outward 
propagating shocks are typically of the order of 100 km/s. 
Some peculiarities arise for the extreme cases of small M u0 
and an infinitely strong incident interplanetary shock (see 
section 3); however, neither of these regimes is significant in 
practice. 
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These results are qualitatively similar to those in Paper 1 
for physically reasonable values of Af w0 , illustrating that the 
conclusions drawn from Paper 1 are robust. It should be 
noted, however, that in Paper 1 an infinite range of p\!p u is 
permitted, whereas in the present paper the corresponding 
range is limited by the jump conditions for the interplanetary 
shocks. The results of the present study suggest that in the 
gasdynamic limit a discontinuous change in pressure can be 
adequately described by a contact discontinuity incident on 
the termination shock. 

In this model, various simplifications have been made. 
One is that the heliospheric magnetic field has been ignored. 
However, as shown in Paper 1, there is an isomorphism 
between the solutions for y = 2 magnetohydrodynamics and 
y =2 gasdynamics, which is valid for any combination of 
adiabatic flows and shocks. Hence the results presented in 
Table 1 for y = 2 would be the same in the MHD case if 17 is 
interpreted as P 3 /P 4 , where P* = p + B 2 !4tt. This result 
is not expected to hold for 7 # 2 or if there is a nonzero 
component of the magnetic field in the flow direction; 
however, the conclusions of the paper are not expected to be 
greatly modified with the inclusion of MHD effects. 

A further simplification has been the neglect of the anom- 
alous cosmic ray component, which Jokipii and Kota [1990] 
suggest could play an important part in governing the struc- 
ture and behavior of the shock. As explained in Paper 1, 
appreciable energy would be used to accelerate the anoma- 
lous component, which would affect the jump conditions and 
shock structure and thickness. Analyses of the even more 
extreme situation in which a shock is modified by galactic 
cosmic rays indicate that the shock structure can be quite 



Figure 3. Plots of V\lv u and V 2 fv u as a function of the 
density ratio p x lp u for 7 = 5/3 and 7 = 2 with M u0 = 100 . 
The velocities are in the rest frame A 0 . V, is the velocity of 
the postinteraction termination shock S x , and V 2 is the 
velocity of the shock S 2 if pj/p M > 1 , or the velocity of the 
head of the rarefraction wave <3l if p\lp u < 1. 



Figure 4. Plots of V\lv u and V 2 lv u as a function of p\/p u 
for various M u0 values with 7 = 5/3. 

broad and complex [e.g., Drury and Volk , 1981; Axford et 
al. y 1982]. Donohue andZank [1993] have recently presented 
a model of the termination shock that incorporates a para- 
metric model of acceleration of the anomalous component; 
their results give an estimate of shock thickness (what they 
call the foreshock) of about 1 AU. It should be noted, 
however, that under the assumptions of the Donohue and 
Zank [1993] model the galactic cosmic rays dominate the 
structure and dynamics of the shock, so that their results do 
not give a clear picture of the situation as it would be if the 
acceleration of the anomalous component should be a large 
effect, while galactic cosmic rays were relatively unimpor- 
tant. 

Another important simplification has been the treatment of 
the shocks as planar, which only gives a reasonable descrip- 
tion of the local and near-term response of the termination 
shock to a change in the upstream dynamic pressure. An 
adequate global description would require generalization of 
the model to at least spherical symmetry. We have also 
assumed that the plasma upstream of the incident shock is 
uniform, whereas in reality it may vary on length scales of 
several astronomical units. A detailed analysis of this situa- 
tion would require numerical simulation. Donohue and Zank 
have simulated such shock pulse collisions, for which the 
postinteraction termination shock evolves back to its origi- 
nal state, for both the gas dynamic and cosmic ray domi- 
nated cases (cf. Figures 8-11 of that paper). As one would 
expect, the region around the postshock contact discontinu- 
ity is considerably more complicated than in the cases 
considered in the present paper. It is not apparent from the 
Donohue and Zank paper whether the interaction results in 
the final termination shock being in motion relative to the 
position of the initial termination shock, as was the case for 
self-canceling density discontinuities discussed in Paper 1 
(cf. Figure 4 of that paper). 
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